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Abstract 

In the framework of a model based on nuclear absorption plus comover inter- 
action, we compute the Et distribution of the J ftp in PbPb collisions at SPS and 
compare it with available NA50 data. Our analysis suggests that the existence 
of new physics (deconfinement phase transition) in the region Et ^ 100 GeV is 
unlikely and that signals of new physics should rather be searched in the region 
Et > 100 GeV. The Et dependence of the J/i/j transverse momentum has been 
computed. At large Et it turns out to be much flatter in the comover approach 
than in a phase transition framework. Estimates of the J ftp suppression at RHIC 
and LHC energies are also given. 
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1 Introduction 



The 1995 data from the NA50 Collaboration show an anomalous J /if} suppres- 
sion, i.e. a suppression larger than the one expected in a nuclear absorption model. This 
model describes the J /if) suppression both in proton-nucleus interactions and in nucleus- 
nucleus interactions with a light projectile Following the original proposal of 0, the 
anomalous J /if) suppression has been interpreted as a signal of a deconfming phase tran- 
sition m |4|, H], H . However, an explanation by a more conventional mechanism, namely, 
the interaction of the J /if) (or the cc pair) with comovers, is also possible |7|, [TD| ] . 



A very spectacular feature of the 1996 data by the same Collaboration |TTJ, is the 
presence of a break in the ratio R(Ej-) of J/if) over Drell-Yan (DY) cross-sections at E? ~ 
55 GeV. It has been argued [|12|, |6| that this break is a signal of deconfmement - although 
there is no general consensus on this point Hl3 , |||, ||, |6j . It is, however, fully recognized 
that a break in the J / ip cross-section would rule out any conventional model, such as the 
one based on comover interaction. (Other approaches to J /if) suppression can be found 

At present, the evidence for this break is weakened by the presence of fluctuations 
in the ratio R(Et) (see Eq. ([T3|) for a precise definition) at large values of Et [ TT[] , which 
are generally regarded as spurious. Also, it is necessary to assess whether this break, if 
confirmed, is due to a genuine break in the J /if) cross-section or rather to fluctuations in 
the J /if) and DY ones. Of course, a definitive answer to these questions can only come 
from data. However, in view of the interest of the subject, it is important to examine the 
available data in a theoretical framework in order to gain some insight on these questions, 
while waiting for a complete analysis of the 1996 NA50 data, and above all, for the 1998 
results. 

The aim of the present work is to perform such an analysis in the framework of a 
model based on nuclear absorption plus comover interaction. This work is a continuation 
of the one in Ref. ||. We use the same formalism and the same values of the parameters 
which were determined in || from the best fit to the ratio R(Et) in pA, SU and PbPb 
collisions. The plan of the paper is the following. In Section 2 we describe the model. In 
Section 3 we compute the Et distributions of minimum bias, DY and J/if>, and compare 
them with available data. In Section 4 we compute the ratio R(Et) and compare it with 
the NA50 data. In Section 5 we compute the Et dependence of the average p T of the J /if) 
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and compare it with recent NA50 data. Section 6 contains our conclusions and prospects 
at higher energies. 



2 The model 



Our model is formulated in a conventional framework |7|, |8], ^, [K| based on two 
different mechanisms of J/ip suppression: nuclear absorption of the pre- resonant cc pair 
with nucleons of the colliding nuclei and absorption by comoving partons or hadrons 
produced in the collision. For completeness we recall its main ingredients: 

Nuclear absorption: In nucleus-nucleus collisions, the survival probability of the J/tp 
at impact parameter b and transverse position s is given by || 

,abs,t. \ t 1 - exp(-A T A (s) a abs )}[l - exp(-B T B (b - s) a ahs )\ 



S abs (b,s) 



a 2 abs AB T A {s) T B (b - s) 

Here TA(B)(b) = f_™ dz PA(B)(b, z) are the nuclear profile functions normalized to unity. 
The nuclear densities pA(B){b, z) are determined from a 3-parameter Fermi distribution 
with parameters given in Ref. fT6| . (In Ref. a different parametrization of the nuclear 
density is used; this introduces differences in R(Et) of less than 4 % in PbPb collisions.) 
For the absorptive cross-section we take a a b s = 6.7 -j- 7.3 mb, consistent with a fit to the 
proton-nucleus data 0. Note that S abs = 1 for a a b s = 0. 

Absorption by comovers. This absorption is due to the interaction of the cc pair (or 
of the J/ip itself) in the dense medium produced in a nucleus-nucleus collision - which 
results in the production of a DD pair. The J/ip survival probability is given by [^. [| 



S c °(b, s) = exp 



/ N co (b s)\ 

-a co N?(b, s) In f 0(Ny°(b, s) - N f ) 



(2) 



Here Ny°(b, s) is the initial density of comovers per unit transverse area d 2 s and per unit 
rapidity at impact parameter b, and Nf is the corresponding freeze-out density. In order 
to have a smooth onset of the comovers, it is natural to take for Nf the density of hadrons 
per unit rapidity in a pp collision, i.e. Nf = [3/(irB%)] dN~/dy\ y * =Q = 1.15 fm" 2 . This 
coincides with the value introduced in Ref. [BJ. With this choice of Nf, the ^-function 
in Eq. (0) is numerically irrelevant. The effect of the comovers in pA turns out to be 
negligibly small. o co is the comover cross-section properly averaged over the momenta of 
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the colliding particles (the relative velocity of the latter is included in its definition). The 
logarithmic factor in Eq. (g) is the result of an integration in the proper time r from the 
initial time to freeze-out time. (One assumes [17, [T3] a decrease of densities with proper 



time in 1/r.) A large contribution to this integral comes from the few first fm/c after the 
collision - where the system is in a pre-hadronic stage. (In this respect, see the last paper 



of Ref. fll0|.) Actually, Brodsky and Mueller [19j introduced the comover interaction 
as a coalescence phenomenon at the partonic level. In view of that, there is no precise 
connection between a co and the physical J/ip — tt or J/ip — N cross-section, and a co has 
to be considered as a free parameter. We take a co = 0.6 mb Q. 

Cross-sections: The J/ip production cross-section in nuclear collisions is given by 

a\ B (b) = [d 2 s m{b, s) S abs (b, s) S co (b, s) , (3) 

Opp J 

where 

m(b,s)=ABa pp T A (s)T B (b-s) . (4) 

We take o pp = 30 mb. With the definition (|3|), the Drell-Yan cross-section (obtained from 
(H) with dabs = Oco = 0) is proportional to AB. 

The cross-section for minimum bias {MB) events is given by 

0AB(b) = 1 - exp[-Op P AB T AB (b)] , (5) 

with T AB (b) =Jd 2 s T A (s)T B (b-s). 

In order to compute these cross-sections we need to know the comover density 
Ny°(b, s) in the NA50 dimuon spectrometer. Moreover, comparison with experiment re- 
quires to compute the above cross-section at a given transverse energy Ej- - measured 
in the NA50 calorimeter. This requires the knowledge of the Et — b correlation function 
P(E T , b). In the following we proceed to calculate these two quantities. 

Density of comovers: It is commonly assumed in the literature that the density of 
comovers is proportional to that of participating (or wounded) nucleons f| ffil . This 
is the so-called Wounded Nucleon Model (WNM; for a review see 0|). In asymmetric 
systems and, in particular, in pA collisions, this model provides a reasonable description 
of the data but only for the average multiplicity - or at negative rapidities, close to the 
maximum of the rapidity distribution. For symmetric AA collisions, the model seems to 
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be valid in a broader rapidity range. (This can be understood from the arguments in 
Ref. see p. 26.) However, for central PbPb collisions (and also for other central 

nucleus-nucleus collisions at SPS) there is experimental evidence of a violation of this 
scaling law at mid-rapidities [^3[ |24|. Moreover, models such as the Dual Parton Model 
(DPM) f22|| , in which unitarity is fully implemented, contain an extra term proportional 
to the average number of collisions. This term is small at present energies but its relative 
size increases with energy. Moreover, it contributes mostly at mid-rapidities. The origin 
of this term is the following. In DPM one has both baryonic strings of type diquark-quark 
and bosonic ones of type q-q. The latter contribute mainly at mid-rapidities. Since the 
number of diquarks available is equal to the number of participating nucleons, the number 
of baryonic strings is equal to the number of participants. On the other hand, the total 
number of strings is proportional to the number of collisions. Therefore, the number of q-q 
strings increases, with increasing centrality, much faster than the number of participants. 
The WNM is obtained from DPM by neglecting the contribution of the q-q strings. 

In the following calculations we will use the density of comovers given by DPM. We 
will also discuss how the J/ip suppression is modified when using a density of comovers 
proportional to the number of participants. 



In DPM, N™(b, s) is given by @, |22 



N™(b, s) = [Ni m A (b, s) + N 2 m B (b, b - s) + N 3 m(b, a)] 6(m B (b, b - s) - m A (b, a)) 
+ [N[ m A (b, s) + N' 2 m B (b, b - s) + m(b,a))6 (m A (b,s) - m B (b, b - a)) . (6) 



Here m is given by Eq. (m and m A , m B are the well known geometric factors |p5 , |XS 



m A{B) (b, s) = A(B) T A(B )(s) 1 - exp (-o- pp B(A) T B(A) {b - s)J . (7) 

The coefficients Ni and N- are obtained in DPM by convoluting momentum distribution 
functions and fragmentation functions [ 2^| . Their values (per unit rapidity) for the ra- 
pidity window and energies of the NA38 and NA50 experiments are given in Table 1 of 
Ref. U . The rapidity density of hadrons is given by 



dN c 



— I d 2 b I d 2 s N™(b,a) 
t arJ J u 



dy a AB 

with <j ab = J d 2 b u AB (b). Note that at fixed b in the range of interest, cr AB (b) ~ 1. 
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The obtained densities of negative hadrons at y* = for pp, SS, SAu and PbPb 
are compared in Table 2 of Ref. || with available data, using in each case the centrality 
criteria (in percentage of total events) given by the experimentalists. In Fig. 1 we compare 
the predictions of both the DPM and the WNM with the NA49 data [|3[] for the rapidity 
distribution of negative particles in central PbPb collisions at 158 AGeV/c. 

Et — b correlation: The experimental results are given as a function of Et- This is 
the total transverse energy of neutrals measured by the NA50 calorimeter in the rapidity 
window —1.8 < y* < —0.6. The correspondence between average values of b and Et is 
given by the proportionality between Et and multiplicity: 

Mb) = q N?Q>) , (9) 

where N™(b) = fl^dy Jd 2 s N™(b,s), with N™(b, s) given by Eq. (|). The parameter 
q is closely connected to the average transverse energy per particle. However, it contains 
extra factors due to the fact that Ny° corresponds to the multiplicity of negatives whereas 
Et is the transverse energy of neutrals. Moreover, a calibration factor of the NA50 
calorimeter (which has an estimated systematic error of about 40 %) is also included in 

A precise determination of q comes from the measured correlation between Et and 
Ezdc ~~ the energy measured at the zero-degree calorimeter. The latter is defined as 

E Z Dc(b) = [A-m A {b)} E in , (10) 

where m^(&) = / d 2 s rriA{b, s), i.e. the average number of participants of A at fixed impact 
parameter, and AE in is the beam energy (E in = 158 GeV/c). A fit to the experimental 
Et — Ezdc correlation using Eqs. (H) and (|10|) allows a precise determination of q. From 
the NA50 data |I| we obtain q = 0.78 GeV. It follows from (|) and ([10]) that with the 



WNM ansatz 0: E T (b) = 0.4 [m A (b) + m B (b)} GeV, the E T - E ZD c correlation is a 
straight line. Experimentally, it is indeed found to be close to a straight line but shows a 
clear concavity. In DPM this correlation has a concavity due to the contribution of the q-q 
strings. However, in the acceptance region of the NA50 Et calorimeter, the contribution 
of the q-q strings is rather small (see Fig. 1) and the concavity is also small. Actually, 
DPM describes well the data in the upper half of the Et region but falls too fast at low 
Et - while the WNM describes the data better in the low Et region (see Fig. 2). One 
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could think that the difference between the two correlation functions is too small to have 
any significant effect on the shape of the distributions. It turns out that this is not the 
case and, therefore, a more accurate description of the Et — Ezdc correlation is needed. 

The failure of the DPM at low Et can be attributed to the effect of the intra- 
nuclear cascade, which is not included in (|6]). This well known phenomenon consists in 
the production of extra particles in the fragmentation regions of the two colliding nuclei 
due to the rescattering of slow secondaries (in the rest frames of the two nuclei) with 
spectator nucleons. Obviously this effect has to vanish for central collisions when no 
spectator nucleons are left. It is also absent at mid-rapidities. However, the rapidity 
region of the Et calorimeter —1.8 < y* < 0.6 is affected by the intra-nuclear cascade 
(which is known to have an extension of about 1.5 rapidity units). In order to incorporate 
the intra-nuclear cascade in a phenomenological way, we replace Eq. (|) by 

E T (b)=qN c y °(b) + kE ZDC (b) . (9') 

With the values of the parameters we use, q = 0.78 GeV and k = 1/4000, the relative 
contribution of the second term in (9') is comparatively small (about 30 % for a very 
peripheral collision with E ZDC = 30000 GeV and less than 2 % for E ZDC < 10000 GeV). 
The only drawback of this extra term is that it does not vanish at E Z dc — ^¥dc = AE in . 
However, this can be easily cured by replacing Eq. (9') by 



E ZDC {b) y E zdc - Ezpcjb) 
4000 / \ E\ 



Mb) = q K W + 0-95 ( ] | ZVC ^max ■ (9") 



J ZDC / 

The corresponding Et — E Z dc correlation, shown in Fig. 2 (full line), is practically 
identical to the one obtained from (9') for Et < 30000 GeV and gives an excellent 
description of the experimental data ||11|| . Moreover, both correlations lead to the same 
Et distributions for J/ip and DY in the region Et > 15 GeV, where data are available. 
Eq. (9") will be used in all DPM calculations. 

In order to obtain the Et — b correlation, and not only the relation between the 
average values of these two quantities, we have to determine the Et distributions at a 
given b. A good description of the experimental E T distributions is obtained |20|, §, 
using a Gaussian distribution at fixed impact parameter, with squared dispersion D 2 {b) = 
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(m°(b)] 2 ) - {N™{b)Y = a(N?(b)), i.e. 



P{E T ,b) 



1 



exp 



[Et - gN y {b)f 
2q 2 aN y (b) 



'27iq 2 aN y (b) 

where N y (b) = E T (b)/q, with E T {b) given by Eq. (9"), and a is a free parameter (see 
Section 3). 



3 Et distributions 

The Et distributions of J/ip, DY and Minimum Bias (MB) are obtained by 
folding the corresponding cross-sections at fixed b (Eqs. (|3]) and (|5|)) with the Et — b 
correlation function: 



dE T 
da MB 



Jd 2 bat B (b)P(E T ,b) , (12) 
= J d 2 ba AB (b) P(E T ,b) . (13) 



dE T 

The corresponding expression for DY is obtained from flT2l) with <r abs = a co = 0. 

The most precise determination of the parameter a is obtained from a fit of (the 



tail of) the MB Et distribution. Using the 1995 data of Ref. [26| we obtain a = 0.73. 
This value will be used in all DPM calculations. With the WNM we use the parameters 
in Ref. f26f : q = 0.4 GeV and a = 1.43. Note that the product aq is the same in both 
cases. (It turns out that the ratio of J/ip over DY is very insensitive to the value of a.) 

The comparison of da DY /dET with the 1995 data |26| is shown in Fig. 3. The 
agreement is satisfactory but the error bars are quite large. Also shown is the distribution 
obtained usine; the WNM. This correlation has a stronger increase with increasing Et ~ 
but both are consistent with the data within errors. 

The comparison with the 1996 Et distribution is shown in Fig. 4. Again the 
(statistical) error bars are quite large. Moreover, there is a significant disagreement both 
with DPM and WNM at Et ~ 135 GeV which was not present when comparing with the 
1995 data. There is also a significant difference in shape between DPM and WNM. Note 
that the only ingredients in the calculation are the b dependence of the DY, ABTAB{b), 
which is common to all models, plus the ET~b or Et — Ezdc correlation. Thus all models 
which reproduce the latter correlation should lead to the same DY distribution. Since 
the DPM (with Eqs. (9') or (9") and (|10"D) gives an excellent description of the latter, 
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the full curve in Fig. 4 should be regarded as the theoretical DY distribution - which 
can be used as a reference when considering the J /if} one. Any significant discrepancy 
with this distribution, such as the one occurring at large Et, should be regarded as a 
possible experimental inconsistency between the measured DY Et distribution and the 
Et — E Z dc correlation. 

We turn next to the Et distribution of the Jfif). We have computed it using 
Cabs — 6.7 mb, a co = 0.6 mb. The result of our calculation is compared with the 1995 



data []26| of the NA50 Collaboration in Fig. 5. The agreement between theory and 
experiment is reasonable. However the data seem to decrease slightly faster than the 
theoretical curve. Note that these data show no break in the J/if> cross-section at any 
value of Et- Fig. 5 shows also the Et distribution obtained with nuclear absorption alone 
(a abs = 7.3 mb) both for DPM (Eq. (9")) and for the WNM. We see that the shape of 
the Et distribution is very sensitive to the effect of the comovers. Thus, a slightly steeper 
decrease of the J/ip cross-section, if confirmed, could possibly be obtained with a small 
increase in the absorption parameters. (The constraint on these parameters coming from 
the SU data is now significantly smaller due to an increase by a factor 2.8 of the statistical 
errors; see Section 6.) Note also that, with nuclear absorption alone, the WNM has a 
faster increase with Et than the DPM. Therefore the extra J /if) suppression required in 
order to reproduce a given shape of the J / if) distribution, must be considerably stronger 
in the WNM than in DPM. This is even more clearly seen in Fig. 6 where we compare the 
theoretical predictions with the J /if) Et distribution from the 1996 data in a linear scale. 
In this comparison we observe some deviations at Et > 100 GeV. This region should be 
studied with great care in the 1998 high statistics run. In our opinion, this is a most 
interesting region to look for eventual signs of new physics, i.e. for the onset of a truly 
anomalous suppression at Et > 100 GeV. On the contrary, in the region Et < 100 GeV 
(where the break in the ratio R(Et) occurs), there is no strong disagreement between 
theory and experiment. However, there is no perfect agreement either and, therefore, it is 
not possible to draw a clear conclusion at present. In particular a sudden drop of the J/if> 



cross-section at Et — 55 GeV has been claimed |TTfl. Even if further data show that this 
drop is statistically significant, our analysis suggests that it cannot be easily attributed to 
a sudden increase of J /if) suppression due to deconfinement. Indeed, in the next four Et 
bins the measured J /if) cross-section is consistent with the predictions of a model which 
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does not have deconfinement. 



4 J ftp over DY ratio 



The J /if) suppression is described by the ratio R(Et) of J /if} and DY cross-sections 
in different Et bins. The advantage of taking this ratio is that systematic errors common 
to both systems do not appear in this ratio. The inconvenient, however, is that the 
results are sensitive to the shape of the DY Et distribution. In our opinion, it is of 
utmost importance to have good data on the Et distribution of the J /if) - as illustrated 
by the analysis of the previous Section. The ratio R(Et) is given by 

K[LT> } d'b a% P(E T ,b) ■ ,14 ' 

This ratio has been calculated, within the present model, in Ref. f9j, where the values 
of the parameters (the same ones used here, including the absolute normalization but 
excepting the value of a which, as discussed in the previous Section, has practically no 
effect on R(Ej-)) were determined from the best fit to R(Et) in pA, SU and PbPb 
collisions. At that time, however, the 1996 data were not available. The comparison of 



the model results with both the 1995 and 1996 data |TjJ is shown in Fig. 7. 

The model reproduces the qualitative behavior of the ratio R. However, there are 
disagreements at a quantitative level. The overall suppression, both from the 1995 and the 
1996 data, is somewhat larger than the theoretical one. More important, the 1996 data 
seem to show a break at Et ~ 55 GeV which is not present in the model calculation. Here, 
several comments are in order. First, as seen in Fig. 7, the experimental data for the first 
Et bin is higher than the one obtained with nuclear absorption alone (with a normalization 
extracted from a fit to pA and SU data [|IJ [|). This is difficult to explain in any model. 
Second, the relevance of the break at Et ~ 55 GeV is weakened by the existence of 
fluctuations in R(Et) at large Et of a comparable size. These fluctuations, which are 
generally regarded as spurious, are an example of systematic errors that do not cancel 
when taking the ratio of J / if) and DY cross-sections and have to be understood. Third, the 
failure of the model to describe quantitatively the ratio R(Et) in the region Et < 100 GeV 
is in sharp contrast with the conclusions reached in Section 3 from a direct comparison of 
the model results with the Et distribution of the J / if) which showed reasonable agreement 



10 



in this Et region. In order to understand the origin of this contradiction we have plotted 
in Fig. 8 the theoretical curve of Fig. 7 (full curve), and compared it with the J/ip 
suppression obtained from the ratio R(Et) of the experimental Et distribution of the 
J/ip (for the 1996 NA50 data) over the theoretical one for the DY (full curve of Fig. 
5). We see that the agreement in shape between theory and experiment has considerably 
improved; the remaining drop, strongly reduced compared with that of R(Et), can be 
attributed to a local minimum in the J/ip cross-section at Et ~ 55 GeV (difficult to 
explain in any model, see comments at the end of the previous Section). The ratio R{Et) 
is now rather well described in the region Et < 100 GeV - except for the first Et bin. 
Moreover, in the ratio R{Et) the break at Et ~ 55 GeV has practically disappeared. Fig. 
8 indicates that the ratio R(Et) is very sensitive to the shape of the DY distribution. We 
would like to stress that Fig. 8 contains no new information. However, it is useful since 
its comparison with Fig. 7 shows the effect on the ratio R(Et) of smoothing out the DY 
cross-section. 

Before concluding this section we would like to comment on the modifications in 
the ratio R(Et) when using the comover density computed in the WNM, rather than the 
one based on DPM, Eq. (|9]). As discussed in Section 2, the WNM underestimates the 
number of negative particles in a central PbPb collision at y* ~ by 15 4- 30 %. The 
corresponding DPM value is 30 % larger than the WNM one and in better agreement with 
the NA49 data (see Section 2). If we would decrease the density of comovers by 30 % for 
the most central Et bin, the value of R(Et) would decrease by about 10 %. Actually, 
the net effect would be significantly smaller, since the WNM multiplicity is smaller than 
the DPM one also in SU, and this can be compensated by a corresponding increase of 
a co . Although a difference would remain, it would not basically change the conclusions of 
the present analysis. On the contrary, our results would be changed if we were to use the 
WNM in the calorimeter rapidity region, in order to determine the Et — Ezdc correlation. 
In this case we would obtain an Et distribution for the J / ip which would be too large in 
the upper half of the Et interval. 
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5 Transverse momentum broadening 



A mechanism producing an increase of the (px) of any type of particle produced in 
pA or AB collisions with increasing nuclear sizes or centrality has been known for a long 
time p7| , |28fl . It is due to initial state rescattering. More precisely, for J/ip production 



this increase is due to rescattering of the projectile and target gluons, before fusion, with 
target and projectile nucleons respectively, encountered in their path through the nuclei. 
The average broadening of the intrinsic gluon distribution in each collision is denoted by 
5q. In Ref. |28| it has been shown that the pt broadening of the J/tp is affected by J/ip 
absorption. In particular the suppression in PbPb collisions in a deconfining approach [[| , 
produces a maximum in the Et dependence of (p T )j/i> & t Et ~ 100 GeV followed by a 
decrease with increasing Et- This peculiar behavior has been considered as a signature 
of Quark-Gluon Plasma formation. 

In this section, we follow the formalism of px broadening of the J / ip in Ref. [ESI , 
but using absorption by comovers instead of the one due to deconfinement. 

The broadening of pr is given by 

5 AB (b) = (p T )A B (b) - (p 2 T ) PP = N AB {b) 5 . (15) 

N AB is the average number of collisions of the projectile and target gluons with target 
and projectile nucleons respectively, up to the formation point of the cc pair, at fixed b. 
This point is specified by the impact parameter b and the positions (s, z) and (b — s, z') 
in the two nuclei. One has 

N AB (b,s,z,z') = a gN A [ dz A Pa{s, z a ) + a gN B [ dz B ps(6 - s, z B ) . (16) 

Here a g N is the gluon-nucleon cross section. This expression (|16|) has to be averaged 
over all positions of the cc formation point with a weight given by the product of nuclear 
densities and survival probabilities: 

W(b, s, z, z') = PA (s, z) p B (b - s, z')S A (s, z) S B (b - s, z') S co (b, s) , (17) 

where 

S A (b, z) = exp (-Aa abs dz p A (b, z)\ (18) 

is the survival probability due to nuclear absorption [15j and S co (b, s) is the survival 
probability due to interaction with comovers, Eq. (EJ). The latter does not depend on the 
cc formation point. 
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We obtain in this way 

_ JdhJ^dzJ^dz' W(b,s } z 7 z>) N AB (b,s,z,z') 
AB[> fd*sf^dz£-dz'W(b,s,z,z>) • 1 ] 

This expression can be written after some transformations as 

J rf 2 s S c °(b, s) [N A (s)D B (b - s ) + N B (b- s)D A (s)] 
NMb) = ° 9N J<Ps S°°(b,s) D A (s) D B (b-s) ' (20) 

where 

D A {s) = —— (1 - exp [-Aa abs T A (s)}) (21) 
Aa abs 

and 

N A (s) = -V (a abs AT A (s) - 1 + exp [-Aa abs T A (s)}) , (22) 

Aa abs 

D B = D A (A -> £) and iV B = A^(A -> 5). 

Finally, the corresponding quantity at fixed transverse energy can be obtained as 

N (F v Jd 2 bP(E T ,b) a AB (b)N AB (b) 
Jd 2 b P{E T) b) a AB (b) 

where P(E T ,b) is the E T — b correlation function, Eq. (|TT|) , and a AB {b) is given by Eq. 
(§• 

The values of (pt)pp an d c^iv^o at 158 AGeV/c are obtained from a fit to the NA50 
data Jll[. One obtains (p^) pp = 1.03^-1.10 (GeV/c) 2 and = 0.39^0.47 ~ 10.0^-12.1 

(GeV fm) 2 (depending on whether the effect of comovers is included or not). The value 
of (t 9N 8q = 0.39 we obtain for nuclear absorption with a abs = 7.3 mb and no comovers, 
agrees with that obtained in [[$|, 9.4 ± 0.7 (GeV fm) 2 , from a fit to pA and SU data 
p^ , [50] . As suggested in PHI , we should take different values of {pj)pp m an d PbPb, 
since this value increases with energy. Using the values measured |52] in 7r~p collisions at 
150 and 200 GeV/c, the value (p%) pp = 1.07 (GeV/c) 2 would correspond to {p^) pp = 1-23 
(GeV/c) 2 at 200 GeV/c. This last value coincides with the one measured in Ref. [32|, 
1.23 ±0.05 (GeV/c) 2 , in pp collisions at 200 GeV/c. (From a fit to DY data in pA, OCu, 
OU and SU [p9] , a value of o^jv^o — 0.13 is obtained, whose ratio over the values for J/tp 
is ~ 0.33, smaller than the value 4/9 ~ 0.44 suggested |27) by the difference of coupling 
between gluons and quarks or gluons.) 

Our results for nuclear absorption plus comovers with a abs = 6.7 mb and a co = 0.6 
mb (nuclear absorption alone with a abs = 7.3 mb), for SU with {p^) PP = 1-23 (GeV/c) 2 
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and a gN 5 = 0.42 (0.40), and for PbPb collisions with {p T )p P — 1-10 (GeV/c) 2 (in agree- 
ment with the mentioned rescaling between 200 and 158 AGeV/c) and the same a g N8o as 
in 577, are shown in Fig. 9 and compared with experimental data JTT], We see that, 
in PbPb collisions with comovers, there is a small maximum at Et ~ 125 GeV. However, 
after this maximum, (p t )ab is practically constant and only slightly smaller than the one 
obtained with nuclear absorption alone. This is in contrast with the sharper decrease at 
large Et found in a deconfming scenario P8 |. The physical origin of this decrease is the 
same in both approaches. At large Et-, corresponding to large comover or energy density, 
the J/ip, with large (p T ) due to a large number of initial g N collisions, is suppressed by 
either the comover or the deconfming mechanisms. However, this effect turns out to be 
numerically much lower in the former approach. Unfortunately, with the present data it 
is not possible to clearly discriminate between these two predictions. 



6 Conclusions and prospects 

We have presented a direct comparison of the available NA50 data for the Et 
distribution of the J/ijj with the results obtained in a conventional framework based on 
nuclear absorption plus comover interaction. 

Our analysis suggests that the presence of new physics in the region Et < 100 
GeV is unlikely. On the contrary, the region Et > 100 GeV is very interesting and should 
be studied with great care in the 1998 high statistics run. Agreement of the J/ip cross- 
section with the comovers results for Et < 100 GeV together with a significantly sharper 
decrease for Et > 100 GeV (for which there might be some hint in the 1996 data), would 
signal the onset of a truly anomalous J/ip suppression. 

Is it possible with the present data to distinguish a deconfming phase transition 
scenario from the more conventional one described here? In order to answer this question 
we have to distinguish between deconfming scenarios producing sharp breaks in the ratio 
R(Et) [12], H from others leading to a smooth behavior of this ratio |3], ^, |5|, |J. For the 



former, a clear-cut answer will probably come from the 1998 data. On the contrary, it will 
be more difficult to distinguish the second type of deconfming models from the comover 
approach presented here. 

A very clear way to do so would be to show that the onset of the anomalous 
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suppression is abrupt, i.e. it is not present below some critical density - for instance the 
maximal one reached in SU collisions ||. Up to a recent date, there was some evidence 



for that [IT]. Indeed, the effect of the comovers in SU produced a somewhat larger 
suppression [RL |9| than the measured one. At present, however, the experimental errors 



in the ratio R(Et) in SU collisions have been increased by a factor 2.8 (see Ref. [[34 



also the experimental errors for PbPb have increased, by a factor 1.4, which has been 
taken into account in this work). In view of that, it is no longer possible to claim that 
the J/ip suppression in SU is too large in the comover approach M or that the onset of 
the anomalous suppression is an abrupt one. 

As we have shown in Section 5, there is a difference between comovers and decon- 
fining scenarios regarding the behavior of (p T ) versus Et- According to Ref. [2^ 



m a 



deconfining scenario this quantity has a maximum at Et ~ 100 GeV and decreases at 
larger Et values. In the comover approach presented here, this drop is practically absent 
and the Et dependence is close to the one obtained with nuclear absorption alone. Al- 
though no such drop is seen in the data, the present experimental errors are rather large 
and a clear conclusion is not possible. 

A promising possibility is the measurement of the J / if) suppression at higher en- 
ergies. The J/ip suppression due to either comover interactions or deconfmement, is 
expected to increase substantially with increasing energy. In the first case, this is due to 
the increase of the density of comovers with increasing energy. In the second case, it is due 
to the corresponding increase of energy density - while the critical value of this quantity 
is unchanged. Therefore, it is important to make predictions at higher energies in both 
approaches, using the values of the parameters determined from present data. One can 
hope that the differences in the predictions of the two approaches will be sufficiently large 
to be experimentally measurable. 

The main uncertainty in the determination of the absolute value of the suppres- 
sion at high energies resides in the value of dN/dy at y* ~ 0. For instance, in central 
PbPb collisions, at RHIC energies, one expects in DPM a value |35] for negative par- 
ticles dN~ J dy\ y *^Q = 1000, and 3500 at ^snn = 5.5 TeV. On the contrary, from the 
scaling in the number of participants (WNM) one expects a value dN~ /dy\ y *^ Q = 400 at 
\J $nn = 200 GeV and 800 at ^/snn = 5.5 TeV. In the first case, there is an increase 
by roughly a factor 5 at RHIC (17.5 at LHC) with respect to the value at y / J/v r /v = 17 
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GeV. In the second case, there is only an increase by a factor 2 at RHIC (4 at LHC), 
which is due to the corresponding increase of dN/dy\ y *^ in pp collisions. An estimate at 
RHIC (LHC) of the J /if} survival probability in central PbPb collisions is given in Table 
1. The numbers in this Table, for comover absorption alone, are obtained from Eq. (0) 
by rising the comover absorption, computed at ^snn = 17 GeV for a central Et bin 
(Et ~ 145 GeV), to a power 2 (4) in the case of the WNM and to a power 5 (17.5) in the 
case of DPM. The corresponding numbers for the total J /if} suppression are obtained by 
multiplying the ones for comovers alone given in Table 1, by the nuclear absorption. The 
latter is expected to depend little on energy [j36fl . 

These estimates illustrate the important increase of the J / if) suppression with ener- 
gy and also the dramatic uncertainties associated to the value of dN/dy\ y *~o. Clearly, a 
more detailed calculation is needed which takes into account the modifications of parton 
densities inside nuclei (usually neglected at SPS energies) and also the changes in the 
Glauber formulae due to the increase with energy of a vv . However, it is obvious that the 
J I if) suppression will increase strongly with increasing energies and it is very unlikely that 
the results will be the same in the comover and in the deconfining frameworks. 
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Table captions: 

Table 1. Comover and total J/ip suppression at SPS, RHIC and LHC, for central PbPb 
collisions, in the WNM and the DPM (see text for an explanation). 
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Figure captions: 



Figure 1. Rapidity distribution of negative hadrons in central PbPb collisions at 158 
AGeV/c. Preliminary data of the NA49 Collaboration (black circles) are compared 
with the DPM results (solid line) using the same centrality criterium, and with the scaling 
in the number of participants (dashed line [f23f). 



Figure 2. — Ezdc correlation: the full line is obtained in DPM from Eqs. (9") and 



TO]) and gives a very good description of the NA50 collaboration data |TT| . The dotted 
line is obtained in DPM from Eqs. (H) and fllQf) . The dashed line is obtained in the WNM 
with E T (b) = 0.4 [m A (b) +m B (b)} GeV §. 



Figure 3. Inclusive cross-section da / dEx for DY pair production with > 4.2 
GeV/c 2 from the 1995 NA50 data [26| compared to the results obtained from Eq. (|TJ) 
with a a iy S = a co = 0. The full curve is obtained in DPM and the dashed one in the WNM. 
The normalization constant a pp Y ja vv in Eq. ([3]) is 9- 10~ 10 . Note that, in order to compare 
with the experimental value of a pp Y /a ppi this normalization factor should be divided by 5 
due to the E? binning in Fig. 3. 



Figure 4. Preliminary E? distribution dN DY / dEx for DY pair production with M w > 
4.2 GeV/c 2 for the 1996 NA50 data |TT[ compared to the theoretical curves of Fig. 3 
normalized to the data. The common normalization factor is 0.10 fm -2 . 



Figure 5. Inclusive cross-section da J ^ '/dE? for J/ip production from the 1995 NA50 
data |2S| compared with the results obtained from Eq. (|T2"|). The normalization constant 
Bw a pl^ I a PP iri © * s ^ ' ^ T ne dotted line is obtained with nuclear absorption 
alone (a a b s = 7.3 mb, a co = 0), while the solid line contains the effect of comovers with 
Cabs = 6.7 mb and a co = 0.6 mb. The dashed line is obtained in the WNM with nuclear 
absorption alone (a a b s = 7.3 mb, a co = 0). 



Figure 6. Premiminary Et distribution dN J ^ /dEx for J/ij) production from the 1996 
NA50 data |TT|, compared with the theoretical curves of Fig. 5 normalized to the data. 
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The common normalization factor is 5.57 fm 2 . The circles and the crosses correspond to 
two different experimental methods [11]: fitting and counting procedures. 



Figure 7. The ratio R(Et) of J/ip over DY versus Et both from the 1995 JIJ (open 
symbols) and preliminary 1996 [[11] (black symbols) NA50 data compared to the ratio of 
theoretical curves (solid lines) in Figs. 4 and 6 (with comovers, solid line). The dotted 
line is obtained in DPM with nuclear absorption alone (cx a fc s = 7.3 mb, a co = 0). The 
normalization factor (61.2) is the one obtained in |9| from a fit to the pA, SU and PbPb 
data. This normalization coincides with the one obtained from the normalizations of the 
individual Et distributions in Figs. 4 and 6 after correcting the latter for the different 
Et binnings, the experimental acceptances and the different DY mass range - which is 
2.9 < M < 4.5 GeV/c 2 in the ratio R(E T ) and M > 4.2 GeV/c 2 in Fig. 4. 

Figure 8. The theoretical curve of Fig. 7 (solid line) is compared to the ratio R(Et) of 
the experimental Et distribution of Fig. 6 over the theoretical DY distribution of Fig. 4 
(solid line). Here the normalization of R{Er) is arbitrary - since we are only interested 
in the change in the shape of R(Et) when smoothing the DY Et distribution. 



Figure 9. (pt)ab in a) SU and b) PbPb collisions at SPS. Solid line: nuclear absorp- 
tion plus comovers; dotted line: nuclear absorption alone (see text for the values of the 
corresponding parameters). Black circles are experimental data from [O, R3|. 
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